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f^ . We propose a new variant model of the modulated reheating. If particles have large scale fluc- 



tuations on their velocities, or equivalently their Lorentz factors, the decay rate also fluctuates 
and the curvature perturbation is induced via their decay processes in analogy with modu- 
lated reheating. For example, if they are produced nonthermally by the decay of another field 
K> , with its mass fluctuating on large scales, such a situation is realized. We explicitly calculate 

C^ I the resulting curvature perturbation and non-linearity parameters and show that the effect of 

velocity-modulation is not negligible if the particles are semi-relativistic at the decay. 



I. INTRODUCTION 

High accuracy measurements of the cosmic microwave background (CMB) anisotropy 
revealed that the cosmological density perturbations obey the nearly scale-invariant power 
spectrum [l|. The inflaton, which drives the inflationary expansion of the very early 
Universe, has quantum fluctuations on large scales and it is a prime candidate for the 
origin of the nearly scale invariant curvature perturbation J2|. 

However, it was recognized that the curvature perturbation can be generated without 
invoking the quantum fluctuation of the inflaton itself. There may be many scalar flelds 
in the physics beyond the standard model, and some of which may be light during in- 
flation and obtain quantum fluctuations. In the curvaton mechanism 3|-[7|, such a light 
scalar other than the inflaton, called curvaton, is responsible for the curvature pertur- 
bation. Another possibility is to make the inflaton decay rate fluctuate due to another 
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light scalar having large scale perturbations [8|, |9[. Then the radiation produced by the 
inflaton decay has curvature perturbation even if the inflaton itself does not have enough 
fluctuations. This is called the modulated reheating scenario. Mechanisms for generating 
the curvature perturbation which shares a similar idea to the modulated reheating have 



been proposed 



m\m. 



In this paper we propose a new mechanism to generate the curvature perturbation, 
which is a variant type of the modulated reheating scenario. The idea is that if a particle 
velocity is fluctuating on large scales, its decay rate also does since the lifetime of a 
particle receives a spatially fluctuating Lorentz boost. Therefore, the modulated reheating 
is realized even if an intrinsic decay rate does not depend on some fluctuating scalar fields 
as in the original modulated reheating scenario. 

As a concrete example, such a large scale fluctuation in the particle velocity is generated 
if it is produced nonthermally by the decay of a heavier particle whose mass depends on 
another light scalar with large scale fluctuations.^ In this setup, modulated reheating 



^ Here the fluctuation of tlie particle velocity should not be confused with the "velocity perturbation" 
usually used in the cosmological context. The former refers to the velocity of each particle while the 



occurs at two stages : one is at the decay of a heavier particle with a fluctuating mass, 
the other is at the decay of a daughter particle whose velocity is fluctuating. The final 
curvature perturbation receives both of the two contributions, as well as that generated 
by the inflaton. We will systematically calculate the resulting curvature perturbation and 
estimate the non-Gaussianity, and find that the effect of modulated velocity is important 
if the velocity at the decay is not much smaller than the speed of light. 

II. BASIC IDEA 

Let us suppose that a particle a has a decay width of T^ at the a rest frame and that 
its velocity is given by v in the laboratory frame. Due to the time delay effect, the effective 
decay rate of a in the laboratory looks like Fg. = T„/'y where 7 = (1 — w^)"^/^. 

Then let us imagine a situation that the a particles are non-interacting and have a 
monochromatic velocity distribution within a small patch of the Universe but on large 
scales their velocities fluctuate. More precisely, we assume that the velocities of cr-particles 
have large scale fluctuations on the decay hypersurface {H{x) = T„{x) in the sudden decay 
approximation where H is the Hubble parameter). On this surface, a has a velocity of 

v{x) = V + 6v{x). (1) 

Finally the cr-particles decay into radiation. Since the decay hypersurface does not coincide 
with the uniform density hypersurface due to the fluctuated Lorentz factor, the curvature 
perturbation is induced at the cr-decay. In this setup, the curvature perturbation is esti- 
mated as 

C ~ nj-^ = -qJ-^ = Q^-f^v6v, (2) 
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where Q„ is the energy density of the a at the decay relative to the total energy density 
and 6r„{x) = Fo-(x) — F^-. Thus the large scale fluctuation on the velocity can be converted 
to the curvature perturbation via the modulated reheating. Obviously, if the a particles 
are non-relativistic {v — ?■ 0), such an effect is negligible. Also in the ultra-relativistic 



latter is the average velocity as the fluid. In the present paper we are interested in the former. 



limit [v — )■ 1), the curvature perturbation is not generated since the equation of state 
of the Universe does not change across the a-decay. Therefore, this effect is expected to 
be important if the a particles are semi-relativistic f ~ (9(1). If the 6v is a Gaussian 
random variable, then the corresponding curvature perturbation is also Guassian at the 
linear order in 6v. 

This is a generic idea and may have a potential of broad applications. In the next 
section we provide a concrete setup which results in the velocity fluctuation and induces 
the curvature perturbation. 

III. A MODEL OF VELOCITY MODULATION AND CURVATURE PERTUR- 
BATION 

A. Outline 

We consider a situation that non-relativistic scalar condensate of the S field decays 
into relativistic a-particles and subsequently cr-particles decay into radiation^. Our key 
assumption is that the mass of S spatially fluctuates. This can be achieved if the mass is 
dependent on a light scalar field which acquired spatial fluctuations during inflation. Such 
fluctuations of the mass result in the generation of the primordial curvature perturbations 
through two different processes. In the following, let us roughly estimate the resulting 
curvature perturbation. More detailed calculations will be given later. 

If the mass of S spatially fluctuates, then generically the decay rate does too. Let us 
suppose that a is a fermion. The S decays into a a pair through the yukawa interaction, 

£ = yT^cra, (3) 

where y is the yukawa coupling constant, then the decay rate is given by 

r(S ^ aa) =T^ = ^ms. (4) 

ovr 



^ The assumption that E is a scalar field is not essential to achieve the velocity modulation. The mech- 
anism can work for the cases where S is not a scalar field as well. 



If (T is a scalar, S can decay into them through the three point interaction, 

L = fiT^aa, (5) 

where /i is a couphng constant having mass dimension one. Then the decay rate is given 
by 

r(E ^ aa) = Ts = -^. (6) 

Svrms 

It is clear from these equations that fluctuation oim-s is taken over by r^. For concreteness 
we assume a is a fermion in the following, although qualitative arguments do not depend 
on whether it is a fermion or scalar. The density fluctuations of a are produced when S 
decays; 
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This is exactly the same mechanism of generating the curvature perturbation in the mod- 
ulated reheating scenario. 

Since each a-particle has an energy m-z/2 when it is produced, velocity of the produced 

cr-particle is given by 

\/ni% — Ami , , 

It is clear that the velocity v also fluctuates if ttie fluctuates. We assume that the interac- 
tion of cr-particle is sufficiently weak so that it does not take part in thermal bath and the 
velocity remains constant except for the redshift by the cosmic expansion.^ Otherwise, the 
cr-particles obey thermal distribution with a background temperature of the radiation and 
the velocity fluctuation will be smoothed out as long as a is subdominant. Now, let us 
suppose that the velocity of a-particle has dropped to vi due to the cosmological redshift 
when (T-particles decay. Taking into account the effect of time delation, the decay rate of 
a is given by 



r((T — )■ radiation) = V „ = V „\ 1 — v1, (9) 



^ The cr-particles can scatter off themselves by the exchange of S. This process is suppressed if the 
coupling constant between S and a is small and/or the S is heavy enough. 



where T„ is the decay rate measured at the rest frame. Since vi is position dependent, 
this equation manifests that the decay rate of a spatially fluctuates too; 

Again, the modulation of To- gives additional contribution to the curvature perturbation. 
Denoting by Q^ a fraction of a-particles when they decay, we expect that the final curvature 
perturbation ( is given by 

c = o> + n,:^ = n,(i + .?)^. (11) 

Pa i^a mj: 

Therefore, if a-particles are still relativistic when they decay, i.e. vi = 0{1), then the 
effect of the velocity modulation on the final curvature perturbation cannot be neglected. 



B. Calculation of the curvature perturbation 



Now let us evaluate the curvature perturbation based on the (5A^- formalism 2l|, |22 |. 
We adopt the so-called sudden decay approximation in which the decay of S or a-particles 
is assumed to occur instantaneously when the decay rate becomes equal to the Hubble 



expansion rate [23 1. 



Let us take a hypersurface just before the S field decays on which the total energy 
density is spatially uniform (uniform density hypersurface). We here assume that E begins 
to oscillate in a period between the end of inflation and the completion of reheating where 
the universe expands like a matter-dominated universe and that the initial amplitude is 
constant. In this case the energy density of the S condensation does not fluctuate on 
the uniform density slice even if its mass fluctuates.^ Since no fluctuations exist before 
S decays, this hypersurface coincides with the spatially flat hypersurface. Let us denote 
by ^A^^i a required e-folding number from the flat hypersurface to the decay hypersurface 



^ If, otherwise, the E begins to osciUate during the radiation-dominated era after the reheating, the 
energy density of E itself also has fluctuations. We do not consider such a case since it only makes the 
following analyses more complicated. 



on which S decays into a-particles. Since F^ spatially fluctuates, 5Ni does too. Thus we 
have 

P.e-^^^^ + -p^e-''^^ = plo, (l + ^) ', (12) 

where pr and ps denote the energy densities of the radiation and S at the S-decay. The to- 
tal energy density is given by ptot = Pr+ps- Quantities with bars represent the background 
values. This is rewritten in the form as 

where ^2^ = pT,/ptot at the S-decay. Let us next consider a hypersurface just after the S 
decays on which po- is spatially uniform and denote by 5N2 a required e-folding number 
from the decay hypersurface to uniform p^. hypersurface. Then equation for 5N2 is given 
by 

where wq is the effective equation of state parameter wq = Pa/ Pa of a evaluated on the 
decay surface. According to the 6N formalism, the sum of 6N1 and 6N2 gives the curvature 
perturbation of the uniform pa hypersurface; 

Ca = SNi + 6N2. (15) 

Eqs. ( IT3l) . ( 1T^ and ( ITSll allow us to express Co- in terms of ST-£, which we defer. Unless 
Wq is zero, (a is non-vanishing, as it should be. Since cr-particles do not interact with the 
other particles, (a is conserved until a-particles decay into radiation. On the other hand, 
(r, the curvature perturbation on the uniform p^ hypersurface, is zero since we assume 
that S decays only into a-particles. 

Now we consider a difference between a hypersurface on which a-particles decay into 
radiation and a one on which pa is uniform. Let us denote by 6N3 a required e-folding 
number from uniform pa hypersurface to a decay hypersurface. Then we obtain 



This is rewritten as 

^^^_3(i+^,)5iV3 ^ (1 _ fi^)e-4(C.+5^3) = A + ^_EA , (17) 

where Q^ = pa/ptot on the a-decay surface and wi is the effective equation of state 
parameter of a on this hypersurface. Let us next consider a hypersurface just after the 
cr-particles decay on which radiation energy density is spatially uniform and denote by 
5Ni a required e-folding number from the a decay hypersurface to the uniform density 
hypersurface. Then we have the following relation 
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g4.yv4 = M + :p j _ (18) 

The final curvature perturbation is given by 

C = C + 5A^3 + 5N,. (19) 

From Eqs. ( TT7|) and ( TTSj) . we can express both SN-^ and 5A^4 as a function of 5Tfj. Therefore, 
C, can be written in terms of 5Ty. and 5Tfj. Note that VLy, and fig- are related through 

f] = ^^^ (20) 

" fis + (l-fis)exp[/3iJ(l + «;(t))rft-4A^]' ^ ^ 

where the e-folding number A^ measures the duration between the S and a decay hyper- 
surfaces and w{t) denotes the time- dependent equation of state of the a particle. The 
integral in the exponent starts from the time of S-decay and ends at the a-decay. It is 
soon seen that in the relativistic limit w = 1/3, Q^ remains constant as is expected. The 
equation of state of a at its production (wq) and decay (wi) are also related. First note 
that the equation of state of the a is given by 

^^') = - = ,r ff^t 2V (21) 

Pa 3(p^(t)+m^) 

because it has a monochromatic momentum distribution in the sudden decay approxima- 
tion, where p{t) is the momentum of the a particle. Then we obtain 



"^1 = n 5 — /sjv I Q — • (22) 

(1 — 3wo)e^ + 3wo 
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Since each a-particle has an energy of ms/2 at the time of their creation, wq can be written 

as 

m% — Ami , , 

To get the final expression for the curvature perturbation, we need to express ^F^ and 
5Y ^ as the functions of ^vhy,- As for the former, it is aheady given by Eq. (j4]). As for the 
latter, the decay rate of a is calculated from 



r,, = r,vi-3wi. (24) 

Through Eqs. fl2^ . fl22|) and fl23l) . we see that Vf, is related to ms- To third order in (5ms, 
we find 



8V a Wi6mY: {wq — 3wi)wi l' 6m-s\ {3wq — 5wi)wl /' 6mj]\' 

Ta Wo rriY. 2wq \ tuy; J 2wl \ mj: J 



(25) 



Using these results obtained above, we can Taylor-expand ( in terms of 6mY. to any 
order. To third order, it is given by 



c=^^+i^.f^y+u(^)', (26) 

ms 2 \mY. J 6 \ mj: J 



where Ai is given by 



^ _ n, [I2wi{wi + 1) - w,{wo + l)(3wi - l)(^s - 4)] 
' 2wo(u'o+l)(fiE-4)((3u;i-l)fi, + 4) ' ^ ^ 

The other expansion coefficients A2 and ^3 are given in the appendix. Let us confirm 

that apart from the 0{1) numerical factors, Eq. ( 1271) reproduces the naive expectation 

Eq. ( ITT1) . First, we see that Ai is proportional to Q^ (It can be confirmed that A2 and ^3 

are also proportional to Q^-)- Therefore, ( of Eq. ( 126|) is proportional to Q^ and vanishes 

in the limit Q^ -^ 0, which is also true for the naive expectation. Secondly, let us expand 

Ai in terms of Wi, 

where we have set Qy. = 0, wq = 1/3 for simplicity.^ We see that the first term remains 
even if we set Wi = 0. Therefore, this term represents the contribution of 6p^ in Eq. 01 ip . 



^ Actually il^r and fls are related through Eq. ([20)1 . Wc can always take fts ^ ^a by choosing N and 

Wo- 



which originates from the modulation of r^. On the other hand, the second term in 
Eq. fl28|) is proportional to Wi and hence proportional to ff. This corresponds to the 
second term in Eq. (ITT]) and is due to the modulation of T„. 



From Eq. f l26|) . the scalar spectral index, n^, is calculated as 24 1 



where we have defined A'^ through Ai5mj]/mY, = N-^5x with x being the light field giving 
spatial modulation to the S mass, Mp denotes the reduced Planck scale, and e and 77 
are inflationary slow-roll parameters [2|. If the infiaton dominantly contributes to the 
curvature perturbation, we recover the standard formula, n^ = 1 — 6e + 2?]. Otherwise, it 
approaches to ra^ = 1 — 2e, similarly to the curvaton case. 

C. Non-linearity parameters 

In this paper, we consider the simplest case where 5ms is a gaussian variable. For 
simplicity, let us assume that S is subdominant when they decay, i.e., VLy, ^ 1 and a- 
particles are relativistic when they are created, i.e., Wq = 1/3. Hereafter we neglect the 
contribution of the infiaton fluctuation to the curvature perturbation and the dominant 
contribution to the curvature perturbation comes from the fluctuation of ms- With these 
assumptions, the non-linearity parameters are given by^ 

20 



/] 



NL 



9(3wi(4w7i - 1) + lyn^iiswi - i)n, + 4) 

X wi{36wi{2wi - 1) + 5) + 1){Q^ - ^WiQ^Y 
+2(wi(3w;i(24w;i(6w;i - 7) + 35) + 10) - 7)n^ 
+2wi{9wi{8wi{Qwi + 5) - 21) + 22) + 22 , 



(30) 



6 



tnl = ^/nl, (31) 

The non-linearity parameters for the standard modulated reheating case are given in Ref. 25[. 
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5'NL 



X 



800 

243(3m;i(4wi - 1) + iyni{-3win^ + ^^ - 4)2 

(3wi - l){wi{3wi{27wi{Awi{96wi{2wi - 3) + 127) - 75) - 43) - 97) + 41)f]^ 
+3{wi{3wi{9wi{wi{12wi{8wi{36wi + 19) - 413) + 2579) - 416) + 262) - 500) + 133)1^2 
+2(u;i(36u;i(3u;i - l)(4u;i - 1) + 1) + l)(fi^ - 3win^)^ 
+2(u;i(3u;i(3w;i(3w;i(12u;i(12(37 - I2wi)wi + 133) - 5177) + 4570) - 728) + 674) - 517)1]<, 



+4(wi + 1)(9m;i(wi(12wi(36wi(4wi + 9) - 181) + 265) - 2) + 209) 

In Fig. [1] and El we show plots of /nl, tnl and (^nl as a function of wi for three cases 
fio- = (0.2, 0.5, 1.0). From this, we see that /nl is always positive. For fixed Wi, the smaller 
Qa is, the larger /nl is. Actually, we can easily show from Eq. f l30|) that /nl — ^a^ when 
Qcr is small. The mechanism of this boost of /nl is exactly the same as that in the curvaton 
model in which /nl is inversely proportional to curvaton fraction evaluated at the time 
when the curvaton decays. For fixed Q^, /nl has a maximum at wi ~ 0.1 and that the 
maximum value is roughly two times larger than /nl at wi = 0. Therefore, if the motion 
of cr-particles are mildly relativistic, then the effect of the velocity modulation can double 

/nl- 

We also find that (^nl is always positive and tnl is always larger than g^i^. Noticing 
that their difference is less than a factor 2 and using Eq. ( 13T1) . we have (^nl = C^(/nl)- 
Therefore, in this model, both tnl and (^nl become large when /nl is large. 



IV. CONCLUSIONS AND DISCUSSION 

We have shown that the curvature perturbation can be sourced by the decay of particles 
with velocity modulation. We have given a concrete setup for realizing such a scenario. 
A class of modulated reheating scenario in which the mass of decaying particle fluctuates 
generates a velocity modulation of daughter particles, which may result in the additional 
source of curvature perturbation. The non-Gaussian signatures may also be enhanced due 
to the velocity-modulation effects. 

One can easily construct a variant model. For example, let us suppose that the mass of 

11 
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FIG. 1: Plots of /nl as a function of wi for three cases r^o- = (0.2, 0.5, 1.0). 

S does not fluctuate but the cr-mass does. The velocity of a, produced by the S-decay, also 
fluctuates in this case and the velocity modulation can generate the curvature perturbation 
at the (T-decay in a same manner. 

The same mechanism may generate the CDM/baryon isocurvature perturbation if the 
cr-decay produces CDM/baryon. Actually a can be identified with the right-handed neu- 
trino Nji while S with a scalar field giving the mass for Nji. The relevant interaction 
Lagrangian is given by 



£ = (y^SiV^^iV^, + yifHL.Nn, + h.c.) - y(S), 



(33) 



where i,j are generation indices and the scalar potential of E is taken to be 



V{^) = -/i||S|2 + A|S| 



(34) 



Here H is the standard model Higgs boson and Lj is the lepton doublet. The S mass 
around the vacuum, ms = 2yUs, is assumed to depend on another light scalar having large 
scale fluctuations. This model possesses a global U(1)l or gauged \J{1)b-l symmetry 
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FIG. 2: Plots of tnl and (^nl as a function of wi for three cases Jla- = (0.2, 0.5, 1.0). 

which is spontaneously broken by the VEV of S.^ In this model S decays into the A'^^;- 
pair, and subsequently Nji decays into the Higgs boson and lepton, generating the lepton 
asymmetry (which is converted to the baryon asymmetry through the sphaleron process) 
if the CP angle is nonzero |27| . Since Nji has a large scale velocity modulation, the baryon 
number created by its decay also has fluctuations. If the E or Nr decay gives dominant 
curvature perturbation, there is no baryonic isocurvature perturbation. On the other 
hand, if the curvature perturbation is dominantly sourced by the inflaton fluctuation, this 
results in the baryonic isocurvature perturbation. 
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Appendix A: Expressions for A2 and ^3 



Here we give expressions for A2 and ^3 in Eq. ([26 



A, 



A, 



Aiw'^{wo + iyi{3wi-i)n^ + Ay 



uwiiwi + lyniiiswi - i)vt, + 4) 

+AAlwlVL„ {wl{wi + 1) {Q?>wlnl - 6wi (7fi(T^ - 2m„ - 8) + 1^1 - 104f]^ + I6O) 
+wl [9wl + Awi + 7) {-?>WiVL„ + f^a - Af + QwoWii'iwi - l){-^WiVL^ + ^a- 4)^ 
+?,wi{?,wi - l)(-3w;ifi^ + ^a- Af) + 'iAiWo^l [Qw^iwi + 1)^ + Uw'lwi {?,wl + 2wi - l) 
+wlwi {27wl + 2Awl + 31u;i - 14) + Qwo{l - ?,wifwl + 3(1 - ?>wifwl) {{?,wi - l)VL^ + 4) 
+wl{'iwi - l)nl [wl (33wJ + 38wi + ^)+wl [9wl + Ibwl + ISwi - 7) + 2wq{1 - ?,Wifwi 

+ (l-3«7i)V) , (Al) 

^,1 + ^,2 + ^,3, (A2) 
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where 



^,1 



3{3wi — l)t/7i(t(7o — 3wi){16AiWo{Qa — 1) + wiQa{—3wiQa + ^o- — 4)) 
"W^o/i wf 3{3wi — l)wlQ^{3wo — 5wi) 



A 



wo + l w^o 2w^o{{3wi - l)n, + 4) 

2 



(A3) 



''' 3wl{wo + l){wi + lyntii^wi - l)fi, + 4) 

-12(^0 + 1)^;;^^ + iYnl{2Aiwo{na - 1) + wiQ^f + 6hwl{n^ - i){win„ + n^y 



l2{wo + l)wl{wi + l)^(fi<, - l)n^^{2AiWo + w,f 

(A4) 



^^'^ = 3^o>o + l)K + lW((3^i-l)l^. + 4) i^"° + '^^"^ + '^^' - "'^)^-(2A,^o + ^i) 
(2Aim;o((3wi - l)fia + 4) + m;i(3wi - l)f],) (2AiWo(wo + l)((3wi - l)fi^ + 4) 
+n^ {w^iwi + 1) + WoWi{3wi - 1) + t/7i(3wi - 1))) {2AiWo{wo + l)((3wi - l)fi^ + 4) 



+ila {Qwl{wi + 1) + woWi{3wi - 1) + wi{3wi - 1))) - 

54f2wl{wo + l){wi + l)5(]^(2Aitf;o(l]^ - 1) + wiQ^) 
{3wi-l)Q^ + 4 

and /i and /2 are defined by 

3u7f (3wi — l)Qa{3wo — 5wi) 



{3wi-l)n^ + 4: 



h 

/2 



2wl{{3wi - l)n^ + A) ' 
3wi(3wi - l)(wo - 3wi)(16AiWo(f^a - 1) + win^{-3win„ + ^^ - 4)) 

2M;3(-3M;if]<, + (^, - 4)2 



(A6) 

■ (A7) 



(A5) 
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